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Abstract —In this paper we analyze and investigate the bit 
error rate (BER) performance of multiple-input multiple-output 
underwater wireless optical communication (MIMO-UWOC) sys¬ 
tems. In addition to exact BER expressions, we also obtain an 
upper bound on the system BER. To effectively estimate the BER 
expressions, we use Gauss-Hermite quadrature formula as well 
as approximation to the sum of log-normal random variables. We 
confirm the accuracy of our analytical expressions by evaluating 
the BER through photon-counting approach. Our simulation 
results show that MIMO technique can mitigate the channel 
turbulence-induced fading and consequently, can partially extend 
the viable communication range, especially for channels with 
stronger turbulence. 

Keywords — MIMO, BER analysis, underwater wireless optical 
communications, log-normal turbulence-induced fading. 

I. Introduction 

Underwater wireless optical communication (UWOC) has 
been recently introduced to meet requirements of high through¬ 
put and large data underwater communications. Acoustic com¬ 
munication systems, which have been investigated and imple¬ 
mented in the past decades have some impediments which 
hamper on their widespread usage for today’s underwater 
communications. In other words, UWOC systems have larger 
bandwidth, lower latency and higher security than acoustic 
communication systems m. These unique features suggest 
UWOC system as a desirable alternative to acoustic commu¬ 
nication systems. 

Despite of all the interesting specifications of UWOC 
systems, they are only suitable for low range underwater 
communications, i.e., typically less than 100 m. This is mainly 
due to the severe absorption, scattering and turbulence effects 
of underwater optical channels. Absorption and scattering of 
photons through propagation under water cause attenuation 
and time spreading of the received optical signals IIl-El. On 
the other hand, underwater optical turbulence which is mainly 
due to the random variations of refractive index (because of 
salinity and temperature fluctuations) results in fading of the 
propagating optical signal a, a. 

Prior works mainly focused on the study of absorption and 
scattering effects of UWOC channels. In IH-a the channel 
turbulence-free impulse response has been simulated and mod¬ 
eled using Monte Carlo simulation method. In 0, a cellular 
topology for UWOC network has been proposed and also the 
uplink and downlink BERs of such a network with optical 
code division multiple access (OCDMA) technique have been 
investigated. On the other hand, some useful studies have been 
accomplished to characterize and investigate turbulence effects 


of UWOC channels. Eor examples, in |0 the scintillation index 
of optical plane and spherical waves propagating in weak 
oceanic turbulence channel, has been evaluated using Rytov 
method. Also the average BER of an UWOC system with 
log-normal fading channel has been investigated in 0, 0- 
Moreover, beneficial application of multi-hop transmission on 
the performance of underwater wireless OCDMA networks has 
been investigated in ifTOll . 

In this paper we analytically study the BER performance 
of an UWOC system, with respect to the all impairing ef¬ 
fects of UWOC channels, namely absorption, scattering and 
turbulence. In order to mitigate turbulence-induced fading and 
therefore to improve the system performance we use spatial di¬ 
versity, i.e., employment of multiple transmitting lasers and/or 
multiple receiving apertures. We assume symbol-by-symbol 
processing and equal gain combining (EGC) at the receiver. In 
addition to evaluating the exact BER, we also evaluate an upper 
bound on the system BER from the inter-symbol interference 
(ISI) viewpoint. Moreover, we use Gauss-Hermite quadrature 
formula and also approximate the sum of log-normal random 
variables with an equivalent random variable, to effectively 
compute the average BERs. 


II. Channel AND System Model 


A. Channel Description 

Propagation of light in underwater medium is under the 
influence of three impairing phenomena, namely absorption, 
scattering and turbulence. Absorption and scattering processes 
cause loss on the received optical signal. Also scattering of 
photons temporally spreads the received optical signals and 
therefore limits the data transmission rate through inducing 
ISI. In order to take into account absorption and scattering 
effects of the underwater channel, we simulate the channel 
impulse response by Monte Carlo simulation method 111-0. 
This turbulence-free impulse response of the UWOC channel 
between any two nodes, (th and jth, is denoted by ho^ij{t). 

On the other hand, turbulence effects of the channel can be 
characterized by a multiplicative fading coefficient, mil- 
ill. Eor weak oceanic turbulence, the aforementioned fading 
coefficient can be modeled as a random variable with log¬ 
normal probability density function (PDE) 0, 0 as; 
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where nxij and . are mean and variance of the Gaussian 

distributed log-amplitude factor Xij = ^\n{hij). Therefore, 
the aggregated impulse response of the channel between any 
ith and jth nodes can be summarized as hij{t) = hijho^ij^t). 
To insure that fading only makes fluctuations on the received 
optical signal, we should normalize fading coefficients as 
E[/iy] = 1, which implies that 

B. System Model 

Consider an UWOC system with M transmitting lasers 
and N receiving apertures. We assume on-off keying (OOK) 
modulation, i.e., the transmitter transmits each bit “1” with 
pulse shape P{t) and is off during transmission of data bit 
“0”. Hence, the total transmitted signal can be defined as 
S{t) = J2 T=-oo bkP{t - kTb), where bk G {0,1} is the fcth 
time slot transmitted data bit and Tb is the bit duration time. In 
the case of transmitter diversity, all the transmitters transmit 
the same data bit bk on their kth time slot. Therefore, the 
transmitted signal of the ith transmitter can be described as 
Siit) = J2’k=-oo^kPt{t - kTb), where = PW^ 

for the sake of fairness. 


A. SISO UWOC Link 

In SISO scheme, the 0th time slot integrated current of the 
receiver output can be expressed a£]; 

-1 
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where h is the channel fading coefficient, = R f^’’ T{t)dt, 
R = ^ is the photodetector’s responsivity, rj is the photode¬ 
tector’s quantum efficiency, q = 1.602 x 10“^® C is electron 
charge, h = 6.626 x 10“®^ J/s is Planck’s constant, / is the 
optical source frequency and L is the channel memory. Fur¬ 
thermore, = R jJ'"’ r(f — kTb)dt = R T(t)dt 

interprets the ISI effect and vt^ is the receiver integrated noise 
component, which has a Gaussian distribution with mean zero 
and variance ifTSl . 

Assuming the availability of CSI, the receiver compares its 
integrated current over each Tb seconds with an appropriate 
threshold, i.e., with Th = /2. Therefore, the conditional 

probability of errors when bits “0” and “1” are transmitted, 
can be obtained respectively as; 


Each ith transmitter, TXi is pointed to one of the receivers. 
The other receivers also capture the transmitted signal of TXi 
due to multiple scattering of photons under water. In other 
words, the transmitted signal of TXi, Si{t) passes through 
channel with impulse response hij{t) to reach the jth receiver, 
RXj. Therefore, the received optical signal from TXi to the 
jth receiver can be determined as; 

00 

= S^{t) * hij{t) = ^ bkhtjrij{t - kTb), (2) 

k——oo 

in which Tij{t) = Piit) * hQ^ij{t) and + denotes convolution 
operation. Furthermore, RXj receives the transmitted signal of 
all the transmitters. Hence, we can express the total received 
optical signal of RXj as; 

M M OQ 
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At the receiver side various noise components, i.e., back¬ 
ground light, dark current, thermal noise and signal-dependent 
shot noise all affect the system operation. Since these com¬ 
ponents are additive and independent of each other, we model 
them as an equivalent noise component with Gaussian distribu¬ 
tion im. We also assume that the signal-dependent shot noise 
is negligible and hence the noise variance is independent of 
the received optical signal (see Appendix A). 

HI. BER Analysis 

In this section we calculate the BER of UWOC system 
for both single-input single-output (SISO) and MIMO con¬ 
figurations. We assume symbol-by-symbol processing at the 
receiver side, which is suboptimal in the presence of ISI 
m. In other words, the receiver integrates its output current 
over each Tb seconds and then compares the result with an 
appropriate threshold to detect the received data bit. In this 
detection process, the availability of channel state information 
(CSI) is assumed for threshold calculation CD. 
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where Q (x) = {1/s/^) ’ exp{—y'^/2)dy is the Gaussian- 

Q function. The final BER can be obtained by averaging the 

conditional BER -p , over 

be\h,bk 2 be\0,h,bk ^ be\l,h,bk 

fading coefficient h and all 2^ possible data sequences for 
bkS, as follows; 
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The form of Eqs. Q and (|6l) suggests an upper bound on 
the system BER, from the ISI point of view. In other words, 
bk^o = 1 maximizes Eq. Q, while Eq. ® has its maximum 
value for = 0- Indeed, when data bit “0” is sent the 

worst effect of ISI occurs when all the surrounding bits are 
“1” (i.e., when bk^o = !)> vice versa lfT4l . Regarding 
to these special sequences, the upper bound on the BER of 
SISO-UWOC system can be evaluated as; 
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^The channel correlation time is on the order of 10“^ to 10“^ seconds 
fs). Therefore, the same fading coefficient is considered for all the consecutive 
bits in Eq. 0 















The averaging in Eqs. 0 and 0 over fading coefficient, 
involves integrals of the form Q{Ch)f^{h)dh, where C is 
a constant, e.g., C = /2(7Tb in second integral of Eq. 0. 

Such integrals can be calculated by Gauss-Hermite quadrature 
formula [15, Eq. (25.4.46)] as follows; 
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in which U is the order of approximation, Wq, q = 1,2, 
are weights of ?7th order approximation and Xq is the gth zero 
of the Uth-order Hermite polynomial, Hu(x) IfTTI . |[T5l . 


B. MIMO UWOC Link 

Assume a multiple-input multiple-output UWOC system 
with equal gain combiner (EGC). The integrated current of 
the receiver output can be expressed as; 


Similar to Section III-A, the upper bound on the BER of 
MIMO-UWOC system can be evaluated by considering the 
transmitted data sequences as bk^Q = 1 for bo = 0 and 
bk^i = 0 for bo = 1. Moreover, similar to Eq. 0 the 
M X iV-dimensional integral in Eq. ([T2l i can be approximated 
by M X W-dimensional series, using Gauss-Hermite quadrature 
formula. 


It’s worth noting that the sum of random variables in 
Eq. O can be effectively approximated by an equivalent 
random variable, using moment matching method Qb). In 
other words, we can reformulate the numerator of Eq. (dUl 
as i-®-’ '^he weighted sum of 

M X N random variables. The weight coefficients are defined 
as = yy + in the special 

case of log-normal distribution for fading coefficients, 
can be approximated with an equivalent log-normal random 
variable as ~ = exp( 22 ;(^“)), with log-amplitude 

mean and variance cr'^^to) of ifTTl : 
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where = R T,^j{t)dt, = R Ti^j{t-kTb)dt = 

R and is the integrated combined 

noise component, which has a Gaussian distribution with mean 
zero and variance A^cr|,^0 

Based on Eq. (ITOl i and availability of CSI, in MIMO 
scheme the receiver selects the threshold value as T/imimo = 
hij^y 12. Therefore, pursuing similar procedures 
as Section III-A results into the following equation for condi¬ 
tional BER. 
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in which H = {hn, hi 2 ,..., Hmn} is the fading coefficients’ 
vector. Assume the maximum channel memory to be Tmax = 
max{Lii, Li 2 , ..., Emm}, then the average BER of MIMO- 
UWOC system can be obtained by averaging over H (through 
M X W-dimensional integral) as well as averaging over all 
2 ^max sequences for &fcs; 
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where /jj(iT) is the joint PDE of fading coefficients in H. 


^Note that the received background power is proportional to the receiver 
aperture area. However based on Appendix A, the background noise has 
negligible contribution on the total noise of the receiver. Moreover, each 
receiver has its distinct dark cun'ent and thermal noise. Hence, the noise 
variance in MIMO scheme is N times of that in SISO case. 


Then averaging over fading coefficients reduces to one¬ 
dimensional integral of 

p(MIMO) f q( \ r Ga^bo)\^ (bo) 

ne\bo,bk -y “ ^ ’ ’ 

which can be effectively calculated using Eq. 0. 

IV. Numerical Results 

In this section we present the numerical results for BER 
of MIMO-UWOC systems. The system is assumed to be 
established in coastal water which has attenuation, absorption 
and extinction coefficients of a = 0.179 m“^, b = 0.219 
m“^ and c = 0.398 m“^, respectively ID. Eurther, we 
assume rectangular pulse shape for transmitted data bits, i.e., 
P{t) = PH j, where P is the total transmitted power 

per bit “1” and n(f) is a rectangular pulse with unit amplitude 
in the interval [—1/2,1/2]. Moreover, we assume the same 
transmitted power of P/M for all the transmitters and the 
same receiving aperture area of A^/N for all the receivers, 
where Ar is the total aperture area of the receiver. Table I 
shows some of the parameters which we have assumed for the 
channel fading-free impulse response simulation and the noise 
characterization (refer to Appendix A for further descriptions 
on the noise components characterization). 

Eig. [T] depicts the upper bound and exact BER of a 25 m 
coastal water link with different configurations and data rate 
of = 1 Gbps. As it can be seen, in all the configurations 
upper bound curves have excellent matches with the exact BER 
curves. Moreover, in a relatively strong turbulent channel, e.g., 
ax = 0.4, spatial diversity (especially at the transmitter side) 
can introduce a noticeable performance improvement, e.g., 6 
dB and 9 dB at the BER of 10“^^, using two and three 
transmitters, respectively. This achievement relatively vanishes 















TABLE I. 


Some of the important parameters for channel 

SIMULATION AND NOISE CHARACTERIZATION. 


Coefficient 

Symbol 

Value 

Half angle field of view 

FOV 

o 

c 

Receiver aperture diameter 

Do 

20 cm 

Source wavelength 

X 

532 nm 

Water refractive index 

n 

1.331 

Source full beam divergence 

Odiv 

0.02" 

Photon weight threshold at the receiver 

■Wth 

o 

1 

a 

Quantum efficiency 

V 

0.8 

Electronic bandwidth 

B 

10 GHz 

Optical filter bandwidth 

AA 

10 nm 

Equivalent temperature 

Te 

290 K 

Load resistance 

Rl 

too a 

Dark current 

Idc 

1.226 X 10“" A 

Received background power 

Pbg 

6.34 X 10“"" W 



Average transmitted power per bit [dBm] 


Fig. 1. Upper bound (UB) and exact BER of a 25 m coastal water link with 
different configurations. = 1 Gbps, cr^ = 0.1 and 0.4. 


in very weak turbulence regimes, e.g., ax = 0.1, where 
fading has a negligible effect on the system performance. 
Furthermore, this figure shows that the transmitter diversity 
performs better than the receiver diversity, due to less noise 
power and larger aperture area. 

In Figs. |2] [2 and ID the upper bound BER of a 25 m coastal 
water link with ax = 0.4 and = 1 Gbps is illustrated 
for different configurations. In particular. Fig. |2] compares the 
results of M X A^-dimensional integrals of Eqs. (fTTI) . (fT2l i 
with the photon-counting method results E). As it is obvious, 
excellent matches between these curves confirm the accuracy 
of the derived expressions in this paper. It is worth noting that 
SIMO schemes, which compensate for fading impairments (at 
high SNR regimes) can not outperform SISO performance, 
except at low BERs. This is mainly due to that each receiver 
in SIMO scheme has N times less aperture area than SISO 
receiver and also SIMO system has N times larger noise 
contribution than SISO scheme. 

In Eig. [2 we applied Gauss-Hermite quadrature formula 
(GHQE) to approximate the M x A^-dimensional integrals 
of Eqs. (fTTI) . (fT^ with M x A^-dimensional series, using 
Eq. @. The order of approximation is assumed to be {7 = 
30. Obviously, GHQE can effectively compute the M x N- 
dimensional integrals (even with less than U = 30 points for 
each integral). In Eig. 0] we used Eqs. (fOT l- lfTSl l to approximate 
the weighted sum of log-normal random variables in Eq. dnii 
with an equivalent log-normal random variable. As it can bee 
seen, this approximation provides an excellent estimate of the 
BER of UWOC system with transmitter diversity. However, 
the discrepancy increases for the case of receiver diversity. 



Fig. 2. Accuracy of the obtained expressions; upper bound BER of a 25 m 
coastal water link with different configurations, /?& = 1 Gbps and crx = 0.4, 
obtained using photon-counting method and classical approach of this paper. 



Average transmitted power per bit [dBm] 

Fig. 3. Comparing the results of M x A^-dimensional series of Gauss-Hermite 
quadrature formula (GHQF) and M x A^-dimensional integrals of Eqs. Hit , 
o in evaluating the upper bound BER of a 25 m coastal water link with 
different configurations, = 1 Gbps and crx = 0.4. 


In other words, in the case of transmitter diversity all the 
transmitters are pointed to a single receiver and therefore, all 
the links have the same weight coefficient of = Gmisq- 
Hence, BER expression of MISO-UWOC can be estimated 
using approximation of unweighted sum of log-normal random 
variables M. 


V. Conclusion 

In this paper we analytically calculated the BER of a 
MIMO-UWOC system with equal gain combining and symbol- 
by-symbol processing. Our analytical treatment included all 
the disturbing effects of the UWOC channels, i.e., absorp¬ 
tion, scattering and fading, we obtained both the exact and 
upper bound BER expressions. Also we used Gauss-Hermite 
quadrature formula to more effectively calculate the averag¬ 
ing integrals with finite series. Moreover, we approximated 
the weighted sum of log-normal random variables with an 
equivalent log-normal random variable to reduce M x N- 
dimensional integrals of averaging (over fading coefficients) to 
one-dimensional integrals. Our analytical results showed well 
match between the exact and upper bound BERs and also the 
results of Gauss-Hermite quadrature formula. Eurthermore, we 
observed that MIMO transmission can introduce a noteworthy 
performance improvement in relatively high turbulent UWOC 



















































Fig. 4. Comparing the results of one-dimensional integral of Eq. (B) and 
M X A^-dimensional integrals of Eqs. E), GD in evaluating the upper bound 
BER of a 25 m coastal water link with different configurations, Rf^ = 1 Gbps 
and (Tx = 0.4. 


channels. However, we assumed log-normal distribution for 
fading statistics, we should emphasize that most of our deriva¬ 
tions are applicable for any fading distribution. 


Appendix A 

Negligibility oe Signal-Dependent Shot Noise 


In this appendix we verify the validity of assumption that 
“the signal-dependent shot noise is negligible with respect to 
the other noise components”. To do that, we should satisfy 
the inequality <C + <^th^ where cr^^, cr|c, 

and are respectively the current variance of the 

Gaussian distributed signal-dependent shot noise, background 
light, dark current and thermal noise QS). In order to verify the 
validity of the above mentioned assumption we should satisfy 
the following inequality; 
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With respect to the parameters in Table I, Eq. dmi simplifies 
to < 6.34 X 10-“ -f 2.688 X 10-®-f 1.097 X 10-3 W Ri 

1 mW. Here, the background noise power is calculated in a 
similar procedure to (dll. Some of the assumed parameters 
are shown in Table I and the other parameters are exactly the 
same as those are in HD. 


To gain more insight on validity of the aforementioned 
assumption let’s to evaluate the BER of an UWOC system with 
pirec) ^ ^ using Gaussian approximation 

ifT^ . Assume an UWOC system with = 10-^ mW 

received power for transmitted bit “1”, which implies to the 

p(rec) 

mean photoelectron counts of nii = —Pfe ~ 2.85 x 10 , 

for Tfc = 1 ns. Also assume the photoelectrons count to be 
Too = rni/2 = 1.425 x 10^, conditioned on transmission of 
bit “0” (mo = mi/2 relates a channel with high ISI). With 
respect to the parameters of Table I, (count) noise variance can 
be obtained as cr^ Ri cr^ rp^j = ~ 3.12 x 10® Il20l . 

Using Gaussian approximation, the BliR of an UWOC system 
with the above parameters for mg, mi and cr^ can be obtained 


as; 
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Therefore, the aforementioned assumption is often valid for 
a wide range of BERs. Note that for typical values of BER, 
p{rec) smaller values than 10-^ W, and hence the above 
assumption is more valid. 
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